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Abstract. We first construct the global unique solution by assuming that the initial data is 
small in the H 3 norm but the higher order derivatives could be large. If further the initial data 
belongs to H~ a (0 < s < 3/2) or (0 < s < 3/2), we obtain the various decay rates of 

the solution and its higher order derivatives. In particular, the decay rates of the density and 
temperature of electron could reach to (1 + in L 2 norm. 



o 

^. \ 1. Introduction 



In the present paper, we consider the compressible non-isentropic Euler-Maxwell system 
(nonconservative form) [JJ [TBI [T6] 

' d t p + div(pu) = 0, 
dtu + u + u- Vu + V9 + GVlnp = -{E + u x B), 

d t Q + u-VQ + ^Gdivu = ^\u\ 2 - (G - 1), 

d t E-V xB = pu, '* ( L1 ) 

d t B + V x E = 0, 
div£ = 1 - p, divB = 0, 
k {p,u,e,E,B)\ t=0 = (p ,u ,e ,E ,B ), x € R 3 . 



The unknown functions p,u,Q,E,B represent the electron density, electron velocity, absolute 
temperature, electric field and magnetic field, respectively. In the motion of the fluid, due to 
■ the greater inertia the ions merely provide a constant charged background. 

Although the compressible Euler-Maxwell system is more and more important in the re- 
l/"") \ searches of plasma physics and semiconductor physics, a small amount of results are obtained 

since its mathematical complexity. In a unipolar form: Chen, Jerome and Wang [2] showed 
the one-dimensional global existence of entropy weak solutions to the initial-boundary value 
problem for arbitrarily large initial data in L°°(R); Guo and Tahvildar-Zadeh |1U| showed a 
blow-up criterion for spherically symmetric Euler-Maxwell system; Recently, there are some 
results on the global existence and the large time behavior of smooth solutions with small per- 
^ ■ turbations, see Tan et al. [23J, Duan [3], Ueda and Kawashima Ueda et al. [27]; For the 

asymptotic limits that derive simplified models starting from the Euler-Maxwell system, we 
refer to |12l [19l [31] for the relaxation limit, |31j for the non-relativistic limit, \17\ [18] for the 
quasi-neutral limit, |24[ [25 j for WKB asymptotics and the references therein. In a bipolar form: 
Duan et al. [4J showed the global existence and time-decay rates of solutions near constant 
steady states with the vanishing electromagnetic field; Xu et al. [32] studied the well-posedness 
in critical Besov spaces. Since the unipolar or bipolar Euler-Maxwell system is a symmetrizable 
hyperbolic system, the Cauchy problem in M 3 has a local unique smooth solution when the 
initial data is smooth, see Kato |14| and Jerome |13] for instance. Besides, we can refer to 
[3 [28] for the non-isentropic case. 
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In this paper, we will refine a global existence of smooth solutions near the constant equilib- 
rium (1, 0, 1, 0, -Bqo) to the compressible non-isentropic Euler-Maxwell system and show some 
various time decay rates of the solution as well as its spatial derivatives of any order. We should 
highlight that our results highly depend on the relaxation terms of velocity and temperature. 
The non-relaxation case is much more difficult, we refer to [6j[8] for such a case. Compared with 
the compressible isentropic Euler-Maxwell system [23J , there are two main difficulties except the 
computational complexity. The first difficulty is that we must obtain the symmetric hyperbolic 
non-isentropic Euler-Maxwell system to do the effective energy estimates. We solve this prob- 
lem by making good use of the positive upper and lower bounds of density and temperature 
(|2.ip . The other difficulty is how to deduce the higher decay rates of density and temperature, 
where is different from the isentropic case [23] since the influence of temperature. To overcome 
this obstacle, we extract two new systems (|3.22p and (|3,28p from the system (|1.2p . then we use 
a bootstrap method to derive (jl.lOp . 

Set 

n = p-l, (9 = 6-1, £ = 5 -.Boo. 
Then the Euler-Maxwell system (jl.ip is reformulated equivalently as 



dtfi = —u ■ Vn — (1 + n)divu, 
d t u + u + E + V0 + u X .Boo = 



Vu 
1, 



1 + 



1 + n 



-Vn — u x B, 



-(l + 0)divu+ 



d t e + e = -u- ve 

dtE — V x B — u = nu, 
d t B |VxB = 0, 
divE 1 = — n, divB = 0, 
_ (n,u,9,E,B)\ t=0 = (n o ,u o ,0o,E o ,B o ), x G 

For N > 3, we define the energy functional by 

TV 

: = || °> E > B ) 

1=0 



(1.2) 



2 

L 2 



and the corresponding dissipation rate by 

N 

V N (t) :=J2\\v l (n,u,0) 



1=0 



N-l 
1=0 



V l E 



L 2 



N-l 

l=i 



\V l B 



L 2 



Our first main result about the global unique solution to the system (|1.2|) is stated as follows. 

Theorem 1.1. Assume the initial data satisfy the compatible conditions 

divEo = —n-o, divBo = 0. 

There exists a sufficiently small 5o > such that if £3(0) < 5q, then there exists a unique global 
solution (n,u,6, E, B)(t) to the Euler-Maxwell system (|1.2|) satisfying 



sup £ 3 (t) + 

0<i<oo 



V 3 (T)dT<C£ 3 (0). 



(1.3) 



Furthermore, if £/v(0) < +00 for any N > 3, there exists an increasing continuous function 
jPjv(') with Pjv(O) = such that the unique solution satisfies 



sup £/v(i) 

0<t<oo 



V N {r)dr < P N (£ N (0)) 



(1.4) 



In the proof of Theorem 1 the major difficulties are the influence of the temperature and 
the regularity-loss of the electromagnetic field. We will do the refined energy estimates stated 
in Lemma l2.8H2.9l which allow us to deduce 

d 



dt 
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—£n + T^n < CnT^n-\£n ■ 
at 



and for N > 4, 



Then Theorem 11.11 follows in the fashion of 1 

Our second main result is on some various decay rates of the solution to the system f) 1 . 2 j) by 
making the much stronger assumption on the initial data. 



Theorem 1.2. Assume that (n,u,9,E,B)(t) is the solution to the Euler-Maxwell system (jl.2p 
constructed in Theorem \l.l\ with N > 5. There exists a sufficiently small So = Sq(N) such that if 
£at(0) < Sq, and assuming that (uo, 9q, Eq, Bo) G H~ s for some s G [0,3/2) or (uo, 0q, Eq, Bq) G 



B 2 ^ for some s G (0, 3/2], then we have 



or 



\\(u,9,E,B)(t)\\^ 3 <C 
\\(u,6,E,B)(t)\\z-s <C . 



Moreover, for any fixed integer k > 0, if N > 2k + 2 + s, then 

V k (n,u,9,E,B)(t) <C7o(l + i)"^- 

L 2 

Furthermore, for any fixed integer k > 0, if N > 2k + 4 + s, then 



V k (n,u,9,E)(t) <C (l + t 



k + l + s 



L 2 



if N >2k + 6 + s, then 
if N > 2k + 12 + s and B, 



L 2 



V k n(t) 
oo — 0, then 
V k {n,9,divu)(t) 



fc + 2 + 3 

<C (l + t) — : 



L 2 



<C Q (l + t) 



(1.5) 
(1.6) 

(1.7) 

(1.8) 
(1.9) 
(1.10) 



In the proof of Theorem 1 1.2 1 we mainly use the regularity interpolation method developed in 
Strain and Guo [21], Guo and Wang |11| and Sohinger and Strain |20| . To prove the optimal 
decay rate of the dissipative equations in the whole space, Guo and Wang [11] developed a 
general energy method of using a family of scaled energy estimates with minimum derivative 
counts and interpolations among them. However, this method can not be applied directly to 
the compressible non-isentropic Euler-Maxwell system which is of regularity-loss. To overcome 
this obstacle caused by the regularity-loss of the electromagnetic field, we deduce from Lemma 
EMU that 



f^_ck+2 
dt k 



+ V k+2 <C k \\(n,u)\\ LOO V k+2 (n,u) V k+ '\E,B) 

L 4 



L 2 



, k aw\x with minimum derivative counts are defined by (|3.9p and (|3.10|) respec- 



where £ k+2 and V k+2 

tively. Then combining the methods of [20] and a trick of Strain and Guo [21J to treat 
the electromagnetic field, we are able to conclude the decay rate (|l,7p . If in view of the whole 
solution, the decay rate (|1.7p can be regarded as be optimal. The higher decay rates (|1.8p - (|1.10p 
follow by revisiting the equations carefully. In particular, we will use a bootstrap argument to 
derive (flTTO]) . 

By Theorem 1 1 . 2 1 and Lemma I2.4H2. 51 we have the following corollary of the usual EP-L? type 
of the decay results: 



2,oo 



Corollary 1.3. Under the assumptions of Theorem \1.2\ except that we replace the H s or B, 
assumption by that (no, 0o, Eo, Bo) G LP for some p G [1, 2], then for any fixed integer k > 0, if 
N > 2k + 2 + s p , then 



V k (n,u,9,E,B)(t) <C (l + t 



fc + Sp 
2 



L 2 



Here the number s 



p •" 



1 
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Furthermore, for any fixed integer k > 0, if N > 2k + 4 + s p> then 



V k (n,u,6,E)(t) <C (l + t) 



p> 

fc+l + sp 



L 2 



(1.11) 



if N > 2k + 6 + s p , then 



if N> 2k + 12 + s p and B 



V k n(t) 
0, then 



L 2 



<c (i + ty 



k + 2 + s„ 



V k {n,6,divu)(t) < C (l + t)-^ + i +s ^. 



L 2 



(1.12) 



The followings are several remarks for Theorem I1.1H1.2I and Corollary 11.31 

Remark 1.4. In Theorem \1. 11 we only assume that the initial data is small in the H 3 norm 
but the higher order derivatives could be large. Notice that in Theorem the H~ s and -B^oo 
norms of the solution are preserved along the time evolution, however, in Corollary \1.3\ it is 
difficult to show that the IP norm of the solution can be ■preserved. Note that the I? decay rate 
of the higher order spatial derivatives of the solution is obtained. Then the general optimal L q 
(2 < q < oo) decay rates of the solution follow by the Sobolev interpolation. 



2,oo 



was introduced there to enhance the 



Remark 1.5. In Theorem \1.2l the space H s or B. 
decay rates. By the usual embedding theorem, we know that for p G (1,2], IP C H~ s with 



1 C B 2 l holds. 



s = 3(- — 2) ^ [0)3/2). Meantime, we note that the endpoint embedding L 1 
Hence the LP-I?{\ < p < 2) type of the optimal decay results follows as a corollary. 

Remark 1.6. We remark that Corollary [L^\ not only provides an alternative approach to derive 
the IP-L? type of the optimal decay results but also improves the previous results of the IP -I? 
approach in Feng et al. [5]. In Feng et al. [5], assuming that B^ = and ||(uo, Eq, Bq)\\ l1 is 
sufficiently small, by combining the energy method and the linear decay analysis, Feng showed 
that 

||(n,0)(i)|| L2 <C (l + f)-^, \\(u,B)(t)\\ L2 <C (l + t)-t and \\E(t)\\ L2 < 0,(1 + t)~i . 

Notice that for p = 1, our decay rate of (n,6)(t) is (1 + t) -13 / 4 in 
in ([LTT]) . 

Notations: In this paper, we use H' 
norm 



fTT2"1) . and u(t) is (l + t) 



-5/4 



s £ R to denote the usual Sobolev spaces with 
and L P (M 3 ), 1 < p < oo to denote the usual LP spaces with norm ||-|| iP . with 
an integer I > stands for the usual any spatial derivatives of order I. When I < or I 
is not a positive integer, stands for defined by tff := J^" -1 (|£|^J^ f), where & is the 



usual Fourier transform operator and ^ is its inverse. We use H s (W 6 ),s G 1 to denote the 
homogeneous Sobolev spaces on M 3 with norm || - 1| defined by ||/||^ 3 := ||A s /|| i2 . We then 
recall the homogeneous Besov spaces. Let 4> £ C^°(i?|) be such that 0(£) = 1 when |£| < 1 
and 0(f) = when |£| > 2. Let tp($) = 0(0 - <p{2£) and tpjfc) = tp(2~^) for j G Z. Then by 
the construction, YljezPjiO = 1 if f 7^ 0. We define Ajf :- 
1 < p, r < 00, we define the homogeneous Besov spaces Bp r 



* /, then for s G R and 
with norm ||-||g s defined by 



Particularly, if r = 00, then 



B s 



LI' 



E 



2 rsj \\&jf\\ r L P 



sup 2 s j Ajf 



LP 



Throughout this paper, we let C denote some positive (generally large) universal constants 
and A denote some positive (generally small) universal constants. They do not depend on either 
k or N; otherwise, we will denote them by Cfc, CV, etc. We will use a < b if a < Cb, and 
a ~ b means that a < b and b < a. We use Co to denote the constants depending on the initial 
data and k,N,s. For simplicity, we write ||(^4,i?)|| x := \\A\\ X + \\B\\ X and J f := L 3 / dx. 
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(*) x e + (**) denote that multiplying (*) by a sufficiently small but fixed factor e and then 
adding it to **. 

The rest of our paper is structured as follows. In section [21 we establish the refined energy 
estimates for the solution and derive the negative Sobolev and Besov estimates. Theorem 11.1 
and Theorem 11.21 are proved in section 

2. Nonlinear energy estimates 

In this section, we will do the a priori estimate by assuming that ||(n, 0)(£)|| H 3 — ^ ^ 
Then by Sobolev's inequality, we have 



- < 1 + n, 1 + < -. 
2 ~ ' ~ 2 



(2.1) 



2.1. Preliminary. In this subsection, we collect the analytic tools used later in the paper. 
Lemma 2.1. Let 2 < p < +oo and a,m,£ > 0. Then we have 

||v a /|| iP <c p ||v m /||i/||vV|[ 2 . 

Here < 9 < 1 (if p = +oo, then we require that < 9 < 1) and a satisfies 

1 I s 



a + 3 



P 



mil 



Proof. For the case 2 < p < +oo, we refer to Lemma A.l in [TT]; for the case p = +oo, we refer 
to Exercise 6.1.2 in [7] (pp. 421). □ 

We recall the following commutator estimate: 

Lemma 2.2. Let k > 1 be an integer and define the commutator 

h = V k {gh) - gV k h. 



V k ,g 



(2.2) 



Then we have 



and 



h 



L 2 



<C k (\\Vg 



v k - l h 



V k (gh) 



L 2 



<c k 



\9\\l° 



V k h 



L 2 



+ 



V k g 



L 2 



\L°° I i 



V k g 



L 2 



Proof. It can be proved by using Lemma 12.11 see Lemma 3.4 in [15] (pp. 98) for instance. □ 

Notice that when using the commutator estimate in this paper, we usually will not consider 
the case that k = since it is trivial. 



Lemma 2.3. If the function f(n) satisfies 



f (k \n) 



< Cfc for any k > 1, 



L 2 



<c k 



L 2 



f(n) ~ n and 
then for any integer k > 0, we have 

\v k f(n) 

Proof. See Lemma 2.2 in |23j. 

We have the IP embeddings: 
Lemma 2.4. Let < s < 3/2, 1 < p < 2 with 1/2 + s/3 = 1/p, then 

\\f\\H-s<\\f\\ LP - 

Proof. It follows from the Hardy-Littlewood-Sobolev theorem, see [7]. 
Lemma 2.5. Let < s < 3/2, 1 < p < 2 with 1/2 + s/3 = 1/p, then 

ll/lk- <II/IL P - 



□ 



□ 
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Proof. See Lemma 4.6 in |20j. 

It is important to use the following special interpolation estimates: 



□ 



Lemma 2.6. Let s > and £ > 0, then we have 



< 



L 2 



v e +lf 



L 2 



where 



1 



l+l + s' 



Proof. It follows directly by the Parseval theorem and Holder's inequality. 



□ 



Lemma 2.7. Let s > and £ > 0, then we have 



V7 

Proof. See Lemma 4.5 in 



< 



L 2 



L 2 



R~ 3 ' 



where 



1 



□ 



2.2. Energy estimates. In this subsection, we will derive the basic energy estimates for the 
solution to the Euler-Maxwell system (|1.2p . We begin with the standard energy estimates. 

Lemma 2.8. For any integer k > 0, we have 

k+2 k+2 



± II V'(n, u, 9,E,B)~ || V'(u, 9) 



dt 



l=k 



2 

L 2 



l=k 



2 
L 2 



'k+2 



a=k 



<C k F(n,u,e,B) [J2\\v\ri,u,0) " 2 +^||V<£ 



fc+i 



l=k 





2 








+ 






L 2 







+ \\(n,u)\\ Loa \\v k+2 (n,u) V k+2 (E,B)\ 

II L 2 I 



L 2 



(2.3) 



where F(n,u,6, B) is defined by 

F(n,u,6,B) := \\(n,u,i 



\ H ^ nH3 + \\(n,u,e)\\ 2 H3 + \\VB\\ L2 . 

Proof. Applying V' (I = k, k + 1, k + 2) to the first five equations in (|1.2p and then multiplying 
the resulting identities by j^V'n, 1 + nV'w, §it^V*0, V'-E, V'i? respectively, summing up 
and integrating over M 3 , we obtain 



1 + 
(it / 1 + n 



| V 'n| 2 + (l + n)|V^| 2 + 31 + n ^^ 12 ' '^^ ^' 2 ' 



2 1 + 



|V'0|' ! + |V l (^ J B)r+ V>,0) 



< 



&0 



1 + 



1 + n (1 + n) z / 
l + f„;, „ „ s „; 31 + n. 



d t n |V z n| 2 + d t n|V'u| 2 + - 



2V1 + (1 + 0) 



L 2 

2 d t e) iv'0^ 



l + n 
1 + 



/ ^-^V z ((l + n)divu)V'n + (1 + n)V ; f i^Vn J • V'u 
J l + n \l + n ) 



v + n v l ev l (\u\ 2 ) 



(i + 0) 



-2 y i±|v , ((l + 0)divw)V'0 + (l + n)V J V0-V , ti + ^ 

-2 J (1 + n)V'(u X J3) • V'u - 2 y nV'tt • V'E - V l {nu) ■ V l E 

:= L 1 +L 2 + L 3 + L 4 + L 5 + L 6 + I 7 . 
First, by (|2.1|) . (jl.2^ . (|1.2|) 3 and Sobolev's embedding inequality, we easily obtain 



(2.4) 



Ji<(l+ \\(n,u,e)\\ H3 )\\(n,u,9)\\ H3 
Next, we estimate the term L 2 . We rewrite L 2 as 



V (n, u, i 



2 

L 2 



I 2 = -2 I v'(u ■ Vn)V l n - 2 [ (1 + n)V'(u • V«) • V'u - 3 /" ^-^V'(u • V0)V z 
J l + n J J 1+9 



hi + I22 + hs- 
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First, we estimate fa- By the commutator notation (|2.2p and (|2.ip . we have 



'21 



l + i 

1 + n 



V'Ut • Vn)V'n 



1 + 



1 + n 



u ■ VV'n + 



< 



u ■ V(VWn) + 



VnV'n 



By integrating by parts, we have 
u ■ V(V'nV'n) 



divu 



Vn 



< ||divu| 



Vn 



L 2 



We employ the commutator estimate of Lemma 12.21 to bound 



V,u 



■ VnV'n < Q ( ||Vit|| LO 



V i_1 Vra 


+ 


Vu 


L2 HVnlUoc) 


Vn 






V 2 








L 2 



<Q||V(n,n)|| io 



V(n,u) 




V z n 






L 2 




L 2 



Then applying the same arguments to fa and /23, by Sobolev's and Cauchy's inequalities, we 
obtain 

7 l, 



L 2 



J2<a||(n,«,e)|| fl3 V'(n,M) 
We now estimate 13. By the commutator notation (j2.2j) . we rewrite ^3 as 



1 + n 

By integrating by parts, we obtain 



2 J( 1 + °) (VdivnV'n + VVn • V z u 

divuV'n + (1 + n) 



1 + 



V',1 + n 



1 + 

1 + n 



Vn • V'w 



2 J( 1 + °) (VdivuV'n + V'Vn • V z -u) = -2 y (1 + 0)div (y l uVn 



2 / V0- VWn< IIV0 



|V0|| L co 






V'n 








L 2 




L 2 



By Lemma 12. 21 and l2.3( Sobolev's and Cauchy's inequalities, we obtain 



1 + ' 



V',1 + n 



1 + n 

<Q||V(n,u,0)|| ioo |v z (n,u, 
In fact, there is a key estimate in the following. 



divnV n + (1 + n) 

2 

L 2 ' 



V 



1 + 

1 + n 



Vn • V z u 



V 



1 + 

1 + n 



Vn 



< C, 



L 2 



V 



1 + 

1 + n 



<Q||V(n,0)|| Lc 

+C z ||(l + 0)|| ioo 
<Ci\\V(n,e)\\ Le 



L°° 



Vn 



V'n 



L 2 



V 



V ; (n, 



L 2 

+ Ci 

1 



V 



1 + 

1 + n 



|Vn| 



L 2 



1 






1 + n 




V ? 







1 + n 



|Vn| 



L 2 



L°° 



<a||V(n,0)|| Zoo V>,0) 



+ Q ||Vn|| Lo 



L 2 



V 1 



IVnl 



1 



1 + n 



L 2 



(2.5) 



where we have used that 1 — ~ n. Then applying the same arguments to I4 and I§, by 
Sobolev's and Cauchy's inequalities, we obtain 

I 3 + h + I 5 < Q || (n,«, 0)11^3 V'(n,n,0) 



L 2 
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For the term Iq, as in |23j . we have that for I = k 



h < C k [ \\u\\ H k + \\VB\\ L2 



V k u 



L 2 



+ 



V k+1 B 



L 2 



(2.6) 



for I = k + 1 



for I = k + 2 



h<C k Midi „ + IIV5 



II 2 



> 




2 

+ 




2 ) 






L 2 




L 2 I 



h ^CfclNI^+^ + llVBlii,, 



) 


V k+2 u 


2 

+ 


V k+1 B 


2 ) 






L' 2 




L 2 1 



+C\\u\ 



V k+2 B 




V k+2 u 






L 2 




L 2 



Next, we estimate the last term Ij. By Lemma 12.21 we easily obtain for I = k or k + 1, 

2 



/ 7 <Q||(n,«) 



Ih 2 







2 




( 


V z (n,u) 


+ 








L 2 





L 2 



for I = k + 2, 



/7<<7 fc ||(n,u)|| ioo V fc+2 (n,n) 



L 2 



v fc+2^ 



L 2 



Consequently, plugging these estimates for l\ ~ I7 into (|2.4[) with Z = fc, k + 1, fc + 2, and 
then summing up, we deduce (|2.3p from (|2.ip . □ 



Note that in Lemma 12.81 we only derive the dissipative estimate of u, 9. We now recover the 
dissipative estimates of n, E and B by constructing some interactive energy functionals in the 
following lemma. 

Lemma 2.9. For any integer k > 0, we have that for any small fixed n > 0, 
'fc+i fc+i 



eft 



^ f V l u ■ VV'n + [ • v ^ - »7 / • V fc V x B 

vZ=fc ^ l=k J 



/k+2 


2 


fe+1 


2 






\z=fc 


L 2 


^\\Ve 

l=k 


+ 

L 2 


V k+1 B 


■-) 



fc+2 



< C^||v'(u,#) ^ + C k G(n,u,6,B) ( ^||v z (n,u,#) 



'fc+2 



vZ = fc 



L 2 



+ 



L 2 



(2.7) 



where G(n,u,9, B) is defined by 



G(n,u,9,B) := \\(n,u,e)f H , +lnH3 + ||VB||£ 2 . 

Proof. We divide the proof into four steps. 
Step 1: Dissipative estimate of n. 

Applying V' (I = k, k + 1) to (|1.2p 2 and then taking the L 2 inner product with VV'n, we 
obtain 



! d t V l u- VVn + /^-^iVV'nl 2 
y 7 1 + n 

< - j 'v l E ■VV'n + C V'(it,V6>) 



1 + n 



L 2 



l + d - Wn 



V l+1 n 



L 2 



+ 



V (u-Vu + ux B) 







L 2 





1 + n 



L 2 

(2.8 
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The delicate first term on the left-hand side of (|2.8p involves dtV l u, and the key idea is to 
integrate by parts in the t- variable and use the continuity equation (jl^p ^ Thus integrating by 
parts for both the t- and x-variables, we obtain 



> 



d_ 

dt 
d_ 
dt 





VV z n 


V l u 


VV'n 


Vu 


VV z n 



V +1 u 



L 2 
2 



J V'divuV' (u ■ Vn + ndivu) 



L 2 



C 



V(u ■ Vn) 



L 2 



c 



V (ndivu) 



L 2 



Using the commutator estimate of Lemma 12.21 we have 



V l (u • Vn) 



L 2 



< 



u ■ V Z Vn 



< u 



L 2 
V m n 



Vn 



L 2 



+ C,\\Vu\ 



<Ci\\(n,u)\\ H3 (\\v l (n,u) 



L 2 



+ 



L 2 

V l+1 n 



+ C, 



L 2 



\Vn\ 



Similarly, 



Hence, we obtain 



( 


V z (n,n) 


+ 


V l+1 u 








L 2 







(2.9) 
(2.10) 



J V l d t u ■ VV l n >j t J V l u • VVn - C 



2 

L 2 



-Ci\\(n,u) 



1^3 



V\n,u) + V l+l (n,u) 



L 2 



L 2 



Next, integrating by parts and using the equation (|1.2p e . we have 



V l E ■ VV'n 



VdivEV l n 



Vn 



L 2 



And as in (|231l ~ (|2TlTl . we have 





+ 


V l+1 u 






L 2 







From the estimate of Iq in Lemma 12.81 we have that for I = k or A; + 1, 

ii 



V (u x B) 



L 2 - V" u H? +1 nH 2 



+ \\VB\\ L2 



)( 


Vu 


+ 










L 2 







(2.11) 

(2.12) 

(2.13) 
(2.14) 



Lastly, by Lemma f2T2l and (|2.5p . we obtain 



V 



1 + 

1 + n 



-Vn 



i±^vv'„ 

1 + n 



V'Vn 



V 



1 



' 1 + n 



Vn- V'Vn<C, || V(n,e)|| Loo ( V\n,9) 



+ 



V m n 



L 2 



.(2.15) 



Plugging these estimates (|2,lip - (|2.15p into ()2.8p . by Cauchy's inequality and (|2,ip . we obtain 

k+i „ fc+2.. .. fc+2 



2 
L 2 



J t E/ V^-VV z n + AE||v z n|[ 2 <C'E||v^,0) 

/fc+2 

+C fc (||(n,M)||^ +lnff3 + ||VB||£ 2 ) ( E \\V l (n, U ,e) 



L 2 



+ 



L 2 



. (2.16) 



This completes the dissipative estimate for n. 
Step 2: Dissipative estimate of E. 
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Applying V' (I = k,k + 1) to (|1.2P a and then taking the L 2 inner product with S/ l E, we 
obtain 



jv l d t u-V l E+ V l E ^ < - J V^|^Vn^ • V l E + C V l (u,V9) 



+ 



1 + n 
V 1 (u-Vu + ux B) 



L 2 



V l E 





V l E 


L 2 





L 2 



L 2 

(2.17) 



Again, the delicate first term on the left-hand side of (|2.17p involves dtV'u, and the key idea 
is to integrate by parts in the t- variable and use the equation (|1.2p 1 in the Maxwell system. 
Thus we obtain 



J V l d t u ■ V l E =j t j W V l E -J VV V l d t E 



d_ 

lit 









V l E- 




W 









V'u • V' (rat + V x B) 



By Lemma 12.21 we have 



V (nu) < Ci\\(n,u)\\ H2 V l (n,u) 



(2.18) 



(2.19) 



We must be much more careful about the remaining term in (|2.18|) since there is no small factor 
in front of it. The key is to use Cauchy's inequality and distinct the cases of I = k and I = k + 1 
due to the weakest dissipative estimate of B. For I = k, we have 



- J V k u ■ V x V k B < e 
for I = k + 1, integrating by parts, we obtain 

-J V fe+ V V x V k+l B = - 

< e 

By Lemma (H2J) 6 , and (|23|) . we have 



V k+1 B 


2 











2 

L 2 



(2.20) 



^ V x V fc+ V V k+1 B 



V k+l B 



L 2 



+ a 



V k+2 u 



L 2 



(2.21) 



i + 



1 + n 



Vn + i^VV'n ) • V z £ 
1 + n / 



v', 1 + 



1 + n 



Vn • V'E 



+ 

< £ 



1 + 9 
l + n 

2 



V z n-V z div£ + I V 



1 + 



l + n 



Vn ■ V l E 



Ve 



L 2 



+ C, £ \\(n,> 



\H3 



V l (n, 



2 

L 2 



(2.22) 



Plugging the estimates (l2T0?]+(|2T2"2l) and (f£I5]) - piljl from Step 1 into (|2TT7|) . by Cauchy's 
inequality, we then obtain 

fc+i „ fc+l .„ „ fc+2 



Z=fc J l=k 



L 2 



< e 



V k+1 B 

k+2 



+ C £ J2 V'(u,fl) 



+CV ( |](n,M)||^ +lnfra + HVSHla) £ |v'(n,n,0) 



;=fc 

2 



2 
L 2 



L 2 



+ 



V k+1 B 



L 2 



. (2.23) 



This completes the dissipative estimate for E. 
Step 3: Dissipative estimate of B . 

Applying V fc to ()1.2p i and then taking the L 2 inner product with —V x \7 k B, we obtain 

2 



V k d t E • V x V k B + 



V x V k B 



L 2 
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< 



X7 k u 



L 2 



V x V k B 



L 2 



+ 



V k {nu) 



L 2 



V x V k B 



L 2 



(2.24) 



Integrating by parts for both the t- and ^-variables and using the equation (|1.2p i . we have 

d 



V k d t E • V x V k B 



dt 
d 
dt 



V k E ■ V x V k B + / V x V k E ■ V k d t B 

2 



- — / V k E • V x V k B 



V x V K E 



L- 



(2.25) 



Plugging the estimates ()2.19p with I = k and f|2.25 j> into f|2.24|> and by Cauchy's inequality, 
since divi? = 0, we then obtain 



- — I V k E-V k V x B + A 
at J 



V k+1 B 



2 

L 2 



< c 


V k u 


2 

+ c 


V k+1 E 


2 






L 2 




L 2 



+ C k \\{n,u)\\ 2 H 2 



V fc (n,u) 



L 2 



(2.26) 



This completes the dissipative estimate for B. 
Step 4 : Conclusion. 

Multiplying (|2.26|) by a small enough but fixed constant rj and then adding it with (|2.23|) so 
that the second term on the right-hand side of (|2.26|) can be absorbed, then choosing e small 
enough so that the first term on the right-hand side of (|2.23|) can be absorbed; we obtain 



d 
dt 



'fc+1 



J V l u ■ V l E -7] J V k E • V fc V x Bj + A ^ || V z 

k+2 2 /k+2 

<CY,\W l (u,0) 2 +C k G(n,u,6,B) ^||V( 



E 



I? 



+ 



V k+1 B 



L 2 



ln,u, 



\l=k 



+ 



V k+1 B 



L 2 



Here G(n,u,9, B) is well-defined. Adding the inequality above with (|2.16p . we get (12 .7p . □ 

2.3. Negative Sobolev estimates. In this subsection, we will derive the evolution of the 
negative Sobolev norms of (u,9,E,B). In order to estimate the nonlinear terms, we need to 
restrict ourselves to that s G (0,3/2). We will establish the following lemma. 

Lemma 2.10. For s € (0,1/2], we have 

| \\(u,e,E,B)f 6 _ s +X\\(u,9)\\ 2 H- s 

(2.27) 



< (J|(n,M)||tf 2 + ||V£||^ ) \\(u,e,E,B)\\^ s + ||£||^ 2 ; 
and for s £ (1/2, 3/2), we have 



d 
~dt 



\\( U ,e,E,B)\\ 2 H _ s + X\\(u,9)\\ 2 H _ 



< I \\(n,u,e)\\ A H i + \\B 



is— 1/2 

\l 2 



VB\\% 2 s \\u\\ L2 ) \\{u,9,E,B)\\ ti . s + \\Ef m . (2.28) 



Proof. The A~ s (s > 0) energy estimate of CLSJ^CLZ]^ yield 

1 d 
2~dt 

9Vn 



(u,E,B)\\%_ a + ^" 2 



A~ s [u ■ Vu + 



1 + n 



+ u x B • A s u 



A~ 



Vn 
1 + n 



■A- S u 



+ 



/a-( 



no) ■ \~*E - I A- s ! ii ■ Y0 + ~9dwu+ l\u\ 2 } A" 



< • Vit + 6>Vn + u x £?| ||it|| + ||nii|| j>_ s ||2£|| jj_ s + || Vnll i, . \\u 



+ ||it • V9 + 9divu + \u\ 



\H-s \\ u \\h- 



3 3' 
I ■ II Ell ■ 



(2.29) 
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We now restrict the value of s in order to estimate the other terms on the right-hand side of 
(l239j) . If s G (0, 1/2], then 1/2 + s/3 < 1 and 3/s > 6. Then applying Lemma EH together 
with Holder's, Sobolev's and Young's inequalities, we obtain 

\W ■ Vu||^_ s < \\u ■ Vu\\ l _l_ j _ < \\u\\ L3/s \\Vu\\ L2 

S MV7 1, 1/2+s |i V7 2 \\ l / 2 - s nv7 II <r II w M 2 . iiT-r ii2 

Similarly, we obtain 

\\9Vn + u x B + nu + u- V0 + 6»divn + |n| 2 ||^_ s < ||(n, u, 0)||# a + ||VS||#i . 

Now if s G (1/2, 3/2), we shall estimate the right-hand side of (|2.29p in a different way. Since 
s G (1/2,3/2), we have that 1/2 + s/3 < 1 and 2 < 3/s < 6. Then applying Lemma |2,4I and 
using (different) Sobolev's inequality, we have 

\\u • Vn||^_ s < ||it|| L 3/ s ||Vn|| i2 < ||«||^2 1/2 l|Vu||^ 2 2_s ||Vu|| L2 

< + ||Vu|| 2 2 ; 

\W x B\\ A _ S < \\B\\ S L2 1/2 \\VB\\% 2 ~ S \\u\y . 

Similarly, we obtain 

\\6Vn + nu + u- V0 + 0divu + |u| 2 ||^_ 3 < \\(n,u,9)\\ 2 Hl . 

Note that we fail to estimate the remaining last term on the right-hand side of (|2.29p as 
above. To overcome this obstacle, the key point is to make full use of the Poisson equation 
(|1.2p 6 . Indeed, using (|1.2p 6 . we have 

l|Vn||£_. < ||A- s Vdiv,E|| L2 < \\E\\ R2 . 

Now collecting all the estimates we have derived, by Cauchy's inequality, we deduce (|2.27[) for 
s G (0,1/2] and {£23} for s G (1/2,3/2). □ 

2.4. Negative Besov estimates. In this subsection, we will derive the evolution of the nega- 
tive Besov norms of (u,6,E,B). The argument is similar to the previous subsection. 

Lemma 2.11. For s G (0,1/2], we have 

±\\(u,9,E,B)\\%. +A||(M)|H. <(\\(n,u,e)\\ 2 H2 + \\VB\\ 2 H1 )\\(u,e,E,B)\\^s + \\E\\ 2 H2 ; 

dZ 2,oo 2,oo\ / A 00 

and for s G (1/2, 3/2] ? we have 

± \\(u,e,E,B)\\%. s +\\\(u,e)\\%- s 

(it D 2,oo D 2,oa 

< (\\(n,u,0)\\w + ll^lll^ 172 HVi?||^ s |M| L2 ) \\(u,6,E,B)\\ KSoo + \\Efjp . 

Proof. The Aj energy estimates of (|l,2P o- (jl.2P c yield, with multiplication of 2~ 2sj and then 
taking the supremum over j G Z, 

5 a(ii^*.^iiki + 5»»iO + l " ll ki + 5»»» , *c. 



< sup 2~ 2sj (- J Aj ■ Vu + Y^Vn + u x B^j ■ Ajuj 

+ 7l 2 ~ 2sj (/ AiM ■ A;/ ' ; - / A ' (irk) • A '") 

+ sup2~ 2 ^' (-|y 4* ( u • v # + ~0divu + ^\u\ 2 ^j ■ AjO^j 



< \\u ■ Vu + 0Vn + u x fill a-s IImIId-s +11^11 d-s II^IId-s +llVn||o- s ||u||o- s 

" n±i 2 ,oo " "-^2,00 " " B 2,ao " n±i 2 ,oo " n±i 2 ,oo " "^2,00 

+ \\u ■ V0 + 0divu + |n| 2 |L- s ||0||r- s . 

11 M± * 2) tx> -"2,00 
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Then the proof is exactly the same as the proof of Lemma 12.101 except that we should apply 
Lemma |2 , 5 1 instead to estimate the B200 norm - Note that we allow s = 3/2. □ 



3. Proof of theorems 

3.1. Proof of Theorem II. 1L In this subsection, we will prove the unique global solution to 
the system (|l,2p . and the key point is that we only assume the H 3 norm of initial data is small. 
Step 1. Global small £3 solution. 

We first close the energy estimates at the H 3 level by assuming a priori that w £ 3 (t) < 5 is 
sufficiently small. Taking k = 0, 1 in (|2.3p of Lemma 12.81 and then summing up, we obtain 



d_ 

dt 



L 2 



1=0 1=0 
Taking k = 0, 1 in (|2.7p of Lemma 12.91 and then summing up, we obtain 



(3.1) 



d 



( E / v ' n ■ vv ' n + E I ylu ■ ylE ~ *i E I ylE ■ v ' v x B J 

\i=o 1=0 J 1=0 J 

/ 3 2 2 \ 3 



2=0 



Z=0 



Since 5 is small, we deduce from (|3.2p xe| ()3.ip that there exists an instant energy functional 
£3 equivalent to £3 such that 

ps + V 3 <0. 

Integrating the inequality above directly in time, we obtain ()1.3p . By a standard continuity 
argument, we then close the a priori estimates if we assume at initial time that £"3(0) < 5o is 
sufficiently small. This concludes the unique global small £3 solution. 
Step 2. Global £n solution. 

We shall prove this by an induction on N > 3. By (|1.3p . then (|1.4p is valid for N = 3. 
Assume (fL4"j) holds for iV - 1 (then now N > 4). Taking fc = 0, . . . , N — 2 in ([23]) of Lemma 
31 and then summing up, we obtain 

N N 



d_ 

dt 



1=0 



1=0 



2 
L 2 



< Cn\/^n-i-\/£n^/T^n + C-^/Pat-i-v/Patv^v < Cn\/T^n-i\/£n\/T^n- (3.3) 

Here we have used the fact that 3 < + 2 < iV — 2 + 1 since N > 4. Note that it is important 
that we have put the two first factors in (|2,3p into the dissipation. 

Taking k = 0, . . . , N — 2 in (|2,7p of Lemma 12.91 and then summing up, we obtain 



d_ 



/N-l r N-l r N-2 . \ 

E / v ' n • vv ' n + E / v ' n ■ y/jE ~ ^ E / ' v x V ' B 

V z=o ^ 1=0 1=0 ^ J 

N-l 



N 



N-l 



+A V V'n + V \\V l E |V V'B 
1 z — ^ II ^2 z — / II r2 z — / 11 



v/=0 

N 



L 2 



1=1 



L 2 



( u > I J L2 + CV y/^N-i V^nV^n- 



(3.4) 



We deduce from (|3.4p x e + ()3.3[) that there exists an instant energy functional £jy equivalent 
to £n such that, by Cauchy's inequality, 

d 



dt 



£n + T^n < CnvDn^iV £nv'E > n < eT>n + Cn^n-i£n- 
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This implies 

— £n + — CnDn-iEn- 

We then use the standard Gronwall lemma and the induction hypothesis to deduce that 

rt t 
S N (t)+ / V n (t)cIt < C£ N (0)e CN fo^- 1 ^^ < C£ N (0)e CNPN - l( - £N - lio)) 
Jo 

< C£ N {0)e CNPN -^ £NW) = P N (€ N (0)) ■ 
This concludes the global £n solution. The proof of Theorem 11.11 is completed. □ 

3.2. Proof of Theorem 11.21 In this subsection, we will prove the various time decay rates of 
the unique global solution to the system (|1.2|) obtained in Theorem 11.11 Fix N > 5. We need 
to assume that £tv(0) < 5o = 5o(N) is small. Then Theorem 11.11 implies that there exists a 
unique global £n solution, and £/v(i) < Pn (£/v(0)) < 5o is small for all time t. Since now our 
5q is relative small with respect to N, we just ignore the ./V dependence of the constants in the 
energy estimates in the previous section. 

Step 1. The H~ s or B^^ norm is preserved along time evolution. 

(|1.5p and (jl.6p indicate that the H~ s or B^^ norm of (u, 9, E, B)(t) is preserved along time 
evolution. First, we prove (|1.5|) by Lemma 12.101 However, we are not able to prove them for 
all s G [0,3/2) at this moment. We must distinct the arguments by the value of s. First, for 
s G (0, 1/2], integrating (pT2Tj) in time, by (fL3|) we obtain that for s G (0, 1/2], 

\\(u,e,E,B)(t)\\l_ s <\\(u ,e ,E ,B )\\%_ s + I P 3 (r)(l+ II («, 9,E,B)(t)\\#-.) dr 

Jo 



<C fl+ sup \\(u,9,E,B)(t)\\^ s ) . 

\ 0<T<t J 



By Cauchy's inequality, this together with (|1.3p gives (jl.5p for s G [0, 1/2] and thus verifies ()1.7p 
for s G [0, 1/2]. Next, we let s G (1/2, 1). Observing that we have (no, 9o, Eo, Bo) £ H^ 1 ^ 2 since 
H- s n L 2 C H- s ' for any s' G [0, s], we then deduce from what we have proved for (|1.7p with 
s = 1/2 that the following decay result holds: 

V k {n,u,e,E,B)(t) <G)(l + i)"^ for A; = 0,1. (3.5) 

H 2 

Here, since we have required N > 5 and now s = 1/2, so we can have taken k = 1 in (|1.7p . 
Thus by (|3.5p . ()1.3p and Holder's inequality, we deduce from (|2.28p that for s G (1/2, 1), 

\\(u,e,E,B)(t)\\%_ s < \\(uo,0o,E o ,B o )\\l_ s + f V 3 (r) (l + \\( U ,9,E,B)(t)\\ A - s ) dr 

Jo 

+ [\\B(t)\\ s - 1/2 \\VB(t)\\% 2 - s V^ \\(u,e,E,B)(r)\\^ a dr 
Jo 

<C fl+fl+ [\l + T)-W- s Mdr) sup \\(u,e,E,B)(r)\\ A . s ) 
V V Jo J o<r<t J 

<C fl+ sup \\{u,e,E,B)(T)\\jj-s). (3.6) 

\ 0<T<t / 

Here we have used the fact s G (1/2, 1) so that the time integral in (|3.6p is finite. This gives 
(fT3j) for s G (1/2,1) and thus verifies ((121) for * G i 1 / 2 , 1 )- Now let s G [1,3/2). We choose 
so such that s — 1/2 < so < 1. Hence, (uo, 0q, Eo, Bq) G H~ Sa . We then deduce from what we 
have proved for (jl.7p with s = so that the following decay result holds: 

V k (n,u,e,E,B)(t) a <C (l + t)- t T SL for k = 0,1. (3.7) 

Here, since we have required N > 5 and now s = so < 1, so we can have taken k = 1 in (|1.7p . 
Thus by (|3.7p and Holder's inequality, we deduce from (|2.28p that for s G [1,3/2), similarly as 
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\\(u,e,E,B)(t)\\%_ 3 <C fl+fl+ f\l + r)-^ 2 -'Ur) sup \\(u,9, E, B)(t)\\ h . 

\ \ J0 J 0<T<t 

<C7 (l + sup \\(u,0,E,B)(t)\U-, J . (3.8) 

V 0<T<i / 

Here we have used the fact s — so < 1/2 so that the time integral in (|3.8p is finite. This gives 
([T5]) for s € [1,3/2) and thus verifies ([177]) for s G [1,3/2). Note that {HD can be proved 
similarly except that we use instead Lemma 12.111 
Step 2. Basic decay. 

For the convenience of presentations, we define a family of energy functionals and the corre- 
sponding dissipation rates with minimum derivative counts as 



k+2 



£ k +2 = Yl \W l (n,u,9,E,B] 



l=k 



2 

L 2 



and 



T>k+2 

k 



k+2 



fe+1 



+ E \\ vlE 



l=k l=k 

By Lemma 12.81 we have that for k = 0, . . . , N — 2, 



L 2 



+ 



V k+1 B 



L 2 



(3.9) 



(3.10) 



d_ 
~dl 



k+2 



E v 



n, tt, 



i=fc 



2 k+2 
),E,B)j L2 +\J2\\V l (u,< 

l=k 
ik+2i 



2 

L 2 



<V^r 2 +||(n, U )|| Lco V fe+ ^(n, U ) V W (£,B) 

! Z, 2 



(3.11) 



By Lemma 12.91 we have that for k = 0, . . . , N — 2, 



'fc+i ~ fc+i ~ » \ 

^ / V'u • VV'n + ^ / V l u ■ V l E - rj V k E • V k V x B 

vZ=fc J l=k J J J 



/k+2 




fc+i 


2 






E 

\l=k 




L 2 ^ II 

;=fc 


+ 

L 2 


V k+1 B 





k+2 2 k+2 



2 

L 2 



(3.12) 



z=fc i=jfc 
Since <5o is small, we deduce from (|3.12p x e + (|3.11|) that there exists an instant energy functional 
£ k+2 equivalent to £% +2 such that 

d 



dt £ k + 2 + V k+2 <\\(n,u)\\ Lac V k+2 (n,u) ^ V k+2 (E,B) 



L 2 



(3.13) 



Note that we can not absorb the right-hand side of (|3.13|) by the dissipation T> k+2 since it does 
not contain ||V fc+2 (-E, B)||^ 2 . We will distinct the arguments by the value of k. If k = or 
k = 1, we bound \\V k+2 (E,B)\\ L2 by the energy. Then we have that for k = 0, 1, 

d 



which implies 



^k+2 + v k+2 < ^k+2^+2^ < ^ +2) 



d 



dt £ k+2 + V k + 2 <0. 



If k > 2, we have to bound || V k+2 (E, B)\\ L2 in term of || X7 k+1 (E, B)\\ L2 since \JV k k +2 can not 
control ||(n,u)|| Loo . The key point is to use the regularity interpolation method developed in 
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[21] . By Lemma \2.1\ we have 



in, U 



V fc+2 (n,u) ( V k+ \E,B) 



L2 



<\\(n,u)\\ L2 2k 



V fc (n,u) 2fc V fe+2 (n,n) : V k+1 (E,B) 



L 2 



where a is defined by 



fc + 2 = (fc + l)x(l 



2A; 



L 2 



+ a x 



L 2 



|V a (£,B)||f 2 ,(3.14) 



2fc 



a = —A; + 1. 
3 



Hence, for k > 2, if N > ffc + 1 <=4> 2<k < |(iV - 1), then by (f3TL4l . we deduce from (13TL3]) 
that 



3/ 



which allow us to arrive at that for any integer A: with < k < '^(N — 1) (note that N — 2 > 



|(iV - 1) > 2 since iV > 5), we have 



d ?fe+2 i ^k+2 
dt k + k 



< 0. 



(3.15) 



We now begin to derive the decay rate from (|3.15|) . Using Lemma [2. 6 1 and (|1.5p . we have that 
for s > and k + s > 0, 





1 




fc + s 




k + s 




L 2 ~ 11 'Iff— 




i2 ^0 


V k+1 B 


fc + l + s 

L 2 



Similarly, using Lemma 12.71 and (|1.6|) . we have that for s > and k + s > 0, 

^ fc-f-s fc-|-5 



V fc B 



< II Rll fc + 1 + s 

J2 ~ 11 "B~ B 
Li 2 , oo 



V fc+i£ 



fc+l+s 

L 2 



< C, 



fc+l+s 

L 2 



On the other hand, for k + 2 < N, we have 

, N-k-2 

-'-!--( / r» ■ < v fc+1 (-E, 5) 
Then we deduce from f)3. 15j) that 

d 



JV_fc - 1 W^iE^B)^- 1 <C V k+1 {E,B) 

L 2 



N-k-2 



N-k-1 
L 2 



^' +2 +K +2 } 1+ "<o. 

where $ = max|^^, N \_ 2 } ■ Solving this inequality directly, we obtain in particular that 



-!/■& 



£ k+2 (t) < { £t +2 (0) +#t> < C (l + t)- 1 ^ = C (l + t)-min{fc+s,JV-fc-2} > (3 16) 



Notice that (|3.16p holds also for k + s = 0ork + 2 = N. So, if we want to obtain the optimal 
decay rate of the whole solution for the spatial derivatives of order k, we only need to assume 
iV large enough (for fixed k and s) so that k + s<N — k — 2. Thus we should require that 



N > max | A; + 2, -k + 1, 2k + 2 + s j = 2k + 2 + s. 

This proves the optimal decay (|1.7p . 
S'tep 5. Further decay. 

We first prove (fL8|) and (fL9j) . First, noticing that -n = divE 1 , by ([LTp . if JV > 2A; + 4 + s, 
then 



V k n(t) 



< 



L 2 



V K+1 E(t) <C (l + t) 

L 2 



(3.17) 



Next, applying V fc to (|1.2|) 2 . p.2|) o. (|1.2|) 1 and then multiplying the resulting identities by 
V fc ii, |V fc #, S7 k E respectively, summing up and integrating over M 3 , we obtain 



2 
1 d 

~2~dt 



V k (u,E) 



2 3 

+ 2 



v k e 



+ 



V fc u 



2 3 
L 2 + 2 



2 

L 2 
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J V k ■ Vu + |^~ Vn + uxBJ ■ V k u + Jv k (V x B + nu) ■ V k E 



- I v k ( u ■ ve + -edwu + -\u\ 2 ) v 



< 



1 + 9 

V k f u ■ Vu H Vn + uxB 



+ 



1 + n 

V fc (n- V6» + 6»divn+ Inl 2 ) 





V fc u 


+ 




L 2 




L 2 





V k (VxB + nu^ 





V k E 


L 2 





L 2 



v k e 



L 2 



L 2 

(3.18) 



On the other hand, taking I = k in (j2. 17j) . we may have 
V fc %/-V fc £+" 



V k E 



L 2 



< 



V k+1 6 


+ 


V k u 


+ 


v fc ( 




L 2 




L 2 





1 + 9 
1 + n 



Vn + it x B 



L' 2 



V k E 



L 2 



.(3.19) 



Substituting (|2.18p with I = k into (|3.19p . we may then have 



d 
~dt 



/ V fc u • V fc £ + 



V k E 



2 

L 2 



< 





2 






+ 


v fc ( 




L 2 





1 + 

1 + n 



Vn laxB 





V k E 


L 2 





L 2 



+ 



v k+l e 


+ 


V k u 


J 


V k E 


+ 




L 2 








L 2 



V k (V x B + nu) 





V k u 


L 2 





L 2 



. (3.20) 



Since e is small, we deduce from (|3.2U|) x e + (|3.18|) that there exists Tk(t) equivalent to 
\\V k (u,6,E)(t)\\ 2 L2 such that, by Cauchy's inequality, Lemma [231 ([23]) . (LTH]) . ([LT]) and ([3+7]) . 

2 



~.F fc (t)+.F fc (t) < 





2 


V k+1 B 


2 




+ 


+ 


v fc ( 




L 2 




L 2 





+ 



V k (u-V9 + 9divu + \u\ 



L 2 



+ 



< 



L 2 



+ (Hull h + ||VB 



+ ||(Vn,VMIl2° 



V k n 



2 

L 2 



V k (nu) 

2 

L 2 



i + e 

1 + n 

2 



Vn + u x B 



L 2 



+ {n,u 



L 2 

V fc+1 £ 



V fc+1 (n,n, 



L 2 



< C7 (l + t) _(fc+1+s) , (3.21) 
where we required iV > 2k + 4 + s. Applying the standard Gronwall lemma to (|3.2ip . we obtain 

-(fc+l+s) 



^(t)<^(0)e-' + C / e-^- T Hl + T)-^ 1+s Ur<C (l + ty 

Jo 



This implies 



V k (u,0,E)(t) r2 <VjW)<C (l + t) 

L 2 



k + l + s 
2 



We thus complete the proof of (|1.8|) . Notice that (|1.9p now follows by (|3.17|) with the improved 
decay rate of E in (|1.8|) . just requiring iV > 2fc + 6 + s. 

Now we prove (jl.lOp . Assuming B^ = 0, then we can extract the following system from 
([I^ - ipreil n. denoting ip = divn, 

5^n + ip = —u ■ Vn — ndivn, 

dtip + i> - n = -A0 - div • Vu + j^ Vn + u x B^ . (3.22) 

Applying V fc to (|3.22j) and then multiplying the resulting identities by V k n, V k ijj, respectively, 
summing up and integrating over R 3 , we obtain 



1 d 
~2dt 



V k n 



+ 



v fc v 



+ 



2 

L 2 



j V k {u ■ Vn + n&\vu)V k n - j V k A9V k 
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j V fc div ■ Vu + Y^ Vn + u x V V ( 3 - 23 ) 



Applying V to (|3.22|) o and then multiplying by — V n, as before integrating by parts over t 
and x variables and using the equation (|3.22p i . we may obtain 



4- I VWn+ 
at 



V fe n 



2 

L 2 

+ I V' 



L 2 



+ / V k nV k ip + I V k (u ■ Vn + ndivu)V k ip 



AO + div [ u • Vu + 



1 + 

1 + n 



Vn luxB 



V fc n. (3.24) 



Since e is small, we deduce from (|3.24p x e + (|3.23p that there exists Gk(t) equivalent to 
V fc (n, VOll 2 j2 such that, by Cauchy's inequality, 



j t Gk(t) + Gk(t) < 



+ 



+ 



V k+1 {u x B) 



V k+1 (u-Vu) 

2 



+ 



L 2 



+ 



L 2 

V k (u- Vn) 



V fe+i 

2 



Vn 



L 2 



+ 



1 + 

1 + n 
V fc (ndivn) 



L 2 



2 

L 2 



(3.25) 



By Lemma 12.21 (|2.5|) and Cauchy's inequality, we obtain 



yfc+1 ( 1 + ^ 
1 + n 



2 

< 






Vn 


2 

+ 


1 + °^k+2 

V k+Z n 


L 2 








L 2 


l + n 



< 



V k+2 n 



L 2 



+ ||V(n,0)||^ V fc+1 (n, 



2 

L 2 
2 

L 2 



and 



V fc+1 (u x B) 



L 2 



u x V fc+1 B + 



2 



V fc+1 ,u 

2 



x B 


2 


« x V fe+1 B 


2 

+ 




V fe+1 ,u 


x B 




L 2 




L 2 









2 

L 2 



L 2 



+ ||V«| 



L°° 



V fc B 



L 2 



+ 



V fc+1 n 



L 2 



IBI 



The other nonlinear terms on the right-hand side of (|3.25|) can be estimated similarly. Hence, 
we deduce from (ET251) that, by (fT77|) - (fTT9]) . 

d :Gk{t)+Gk{t) 



df 



< 



V k+2 (n,0) 

2 



L 2 



+ u 



V fc+1 B 

2 



L 2 



+ ||V«| 



V fc B 



L 2 



+ \\B\\ 



V fc+1 u 



L 2 



V k+2 (n,u) +\\V(n,u,0)\\ 2 Loo V k+1 (n,u,9) 

L 2 



L 2 



+ IK™,«)lli/*> 

< Co ((1 + t)-( fe+3+s ) + (1 + t)-(*+7/2+2 8 ) + (1 + t y(k + n/2+2s) j 

< C (l + t)" (fc+3+s) , (3.26) 
where we required iV > 2k + 8 + s. Applying the Gronwall lemma to (|3.26j) again, we obtain 



G k {t) <&(0)e-* + <7 / e 

JO 



-(*- T )(l + T )-( fc+3+s ) dr < C (l + t)-( fc + 3 + s ). 



This implies 



fc+3+s 



V*(n,^)(t) r < 0,(1 + t) 2 

L 2 



(3.27) 



Now we consider the following system which consists of (|3.22p and (|1.2P o: 

dtn + if) = —u ■ Vn — ndivu, 

(1 + 
u ■ Vu H Vn + n x B 
l + n 

d 4 + = -u • V0 - \(l + 0)divu + \\u\ 2 . 



(3.28) 
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First, we have the standard energy identity for the system (|3.28p 



1 d 

2 ~dt 



V k n 



+ 



v fe v 



+ 



+ 



L 2 



J V fc (u • Vn + nd\wu)V k n - j V fe A0VV - | J V k ^V k t 
V fe div ( u ■ Vu + + u x B J V fc V 



1 + n 



V fe u • V6» + -0divu 



d 2 V fe 6 



(3.29) 



Since e is small, we deduce from (|3.24p x e + (|3.29p that there exists % k {t) equivalent to 



| V k (n, tp, 6) |L 2 such that, by Cauchy's inequality and ()3.27p 



-H k (t) + H k (t) < 



L' 2 



+ 



V fc+1 (u • Vu) 



L- 



V fc+1 (u x 5) 



+ 



v fc v 



L 2 



+ 



+ 

L 2 

fc/i„.i2 



V k (u ■ Vn. 



+ 



L 2 



+ 



L 2 

V k (u-V0) 



1 + 71 

V fe (radivu) 

2 



2 

L- 



2 

L 2 



L 2 



+ 



V fc (0divu) 



L 2 



v*(i«r) 

< C ((1 + t)~(*+3+«) + (1 + t )-(fc+7/2+2s) + (1 + ^-(fc+11/2+2^ 

< C (l + t)~ (fc+3+s) , (3.30) 
where we required N > 2/c + 8 + s. Applying the Gronwall lemma to (|3,30p again, we obtain 



n k {t) < n k {o)e-* + c 

This implies 



<t-r) 



(l + r)-( fc+3+s ) dT<C (l + t) 



-(fc+3^ 



V k (n,6,m) <V^)<Co(l + *) 

L 2 



(3.31) 



If required N > 2k + 12 + s, then by ()3.31 j) . we have 

7 k+2 



V k+ \n,9){t) <0,(1 + 

L 2 

Having obtained such faster decay, we can then improve ()3.26p to be 

j t Q k {t) + G k (t) < C ((1 + t)-( fc +5 +S ) + (1 + tr (k+7/2+2 S ^ < Co(1 + i) _( fc+7 /2+2.) i 



Applying the Gronwall lemma again, we obtain 

V k (n,4,)(t) < 



L 2 



VGhV) < 0,(1 + 1) 



-(fc/2+7/4+s) 



In light of the faster decay for V fc V) we can then improve ()3.30p to be 

jU k {t) + U k {t) < Co ((1 + t)-( k+5+s ^ + (1 + t )-(k+7/2+2s)^ < + 
Applying the Gronwall lemma again, we obtain 

V k (n,9,m) a < < C (l + t)^( fc / 2 + 7 /4+ S ). 

L 2 

We thus complete the proof of (jl.lQj) . The proof of Theorem 11.21 is completed. 
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